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Hints of direct dark matter detection coming from the DAM A, CoGeNT 
experiments point toward light dark matter with isospin-violating and pos- 
sibly inelastic couplings. However an array of astrophysical constraints are 
rapidly closing the window on light dark matter. We point out that if the 
relic density is determined by annihilation into invisible states, these con- 
straints can be evaded. As an example we present a model of quasi-Dirac 
dark matter, interacting via two U(l) gauge bosons, one of which couples 
to baryon number and the other which kinetically mixes with the photon. 
Annihilation is primarily into "dark neutrinos" that do not mix with the SM, 
but which could provide an extra component of dark radiation. The model 
could soon be tested by several experiments searching for such light gauge 
bosons, and we predict that both could be detected. The model also requires 
a fourth generation of quarks, whose existence might increase the production 
cross section of Higgs bosons at the Tevatron and LHC. 

Introduction. The DAMA [Ij and CoGeNT ^ Q\ experiments have presented evi- 
dence for light, ~ 10 GeV dark matter (DM), which is at odds with null results from 
XenonlO 0], XenonlOO [5] and CDMS [6j for the simplest DM models, and moreover the 
cross sections needed by DAMA and CoGeNT are at odds with each other. It is very 
intriguing that a single hypothesis, that DM has isospin-violating interactions with nu- 
cleons, resolves the discrepancies between all of the experiments except CDMS [7j-[12j. 
One further assumption, that the interactions are inelastic, connecting two DM states 
split by ~ 10 keV, helps to alleviate the tension with CDMS [HIS] (see also [I3])Q 
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^The CRESST experiment [13] reports evidence for light DM with mass 9 GeV at the lower end of 
their M2 best fit region and cross section an ~ 1.5 x 10^* pb. In an isospin violating model with 
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Recently we proposed some "minimal models" of hidden sector dark (DM) matter 
[15] that have the desired properties, and we noted that a number of astrophysical 
considerations constrain the models very strongly. Essentially, if DM of mass ~ 10 GeV 
annihilates primarily into any channel other than muons, it is ruled out by constraints 
from dumping electromagnetic energy into the cosmic microwave background (CMB) 
[Ml E], Fermi observations of dwarf satellite galaxies [Tgl [TS], or SuperK limits on 
neutrinos from DM annihilations in the sun |20^ [2T| [22] . These constraints will tighten as 
a result of forthcoming data from new experiments like Planck. Moreover the PAMELA 
constraint on cosmic ray antiprotons excludes 10 GeV DM with an annihilation cross 
section greater than 0.1 times the standard relic density value if the final state contains 
quarks that can hadronize |23j . This tension is demonstrated to be rather insensitive to 
different choices for cosmic ray propagation models and halo models for the bb channel 
in [24]. The only robust particle physics mechanism for evading this tension is if the 
annihilation is into a pair of bosons that are too light to decay into p-p (or do not couple 
to quarks). 

Of course one of these new experiments may find deviations giving positive indirect 
evidence for light dark matter. But if they instead only tighten the constraints, while 
direct evidence for light DM persists, we will have a puzzle. One elegant way out is 
asymmetric dark matter that carries a conserved charge. If the symmetric component 
of the DM has annihilated away, then the above constraints no longer apply. In the 
present letter, we wish to point out a different possibility, assuming the DM is symmetric. 
Namely, if it annihilates primarily into invisible particles, the constraints in question are 
evaded. While this may seem like a trivial statement, in fact care must be taken to 
avoid problems from other constraints, such as the number of species during big bang 
nucleosynthesis (BEN) or too-rapid cooling of supernovae by emission of the new invisible 
statesn 




To make these points concrete we will illustrate them in a hidden sector model similar 
to one presented in [R]- This model achieves isospin violation using a light vector 
boson that couples both to gauged baryon number and, and another one Z' that 
kinetically mixes with the photon. The DM gets inelastic couplings and a small mass 
splitting through a weak Yukawa interaction with a new Higgs field that breaks the U(l) 
gauge symmetries. It annihilates dominantly into light "dark neutrinos" via exchange 
of another new Higgs field, which is also the one primarily responsible for the ~ 10 GeV 
DM mass. Thus the model is economical, in that most of elements serve for more than 
one purpose. An interesting feature of this model is that the kinetic mixing of the Z' 
boson is related to its mass in such a way that it might be discovered in the near future 
by several proposed experiments aimed at detecting such particles. Moreover because 
the B boson might mix weakly with Z' , it could also be found in the same searches. 



fp/fn = —1.5, this value would increase by a factor of + fp/UT^ = 16, since O and Ca have 
equal numbers of protons and neutrons, which is a factor of 13 below the inelastic CoGeNT/DAMA 
best-fit value for o"„ = 0.03 pb which we adopt below [11]. This discrepancy is greatly reduced (to 
a factor of 2) if one adopts the value an — 5 x 10~^ pb, corresponding to the best fit with elastic 
scattering. 

^We thank D. Spolyar for reminding us about this important constraint. 




2 



The Model. Gauged baryon (B) number is anomalous in the standard model, but the 
anomalies can be canceled by adding a fourth generation of quarks with baryon number 
±1 [25l[26]. If we wish to further couple the DM x to the associated vector boson S^, 
then to avoid further anomalies it is natural to assume that x is initially vector-like 
(Dirac), but a new Higgs boson (j) that spontaneously breaks some combination of the 
dark U(l) symmetrires can also couple to x ^^id render it quasi-Dirac after symmetry 
breaking. This fits in nicely with the additional preference that the DM interactions 
with nuclei should be inelastic. A small ~ 10 keV mass splitting between the two DM 
components is thus both welcome and natural. But because the preference for inelasticity 
is weak in global fits to all data [TT] . we will consider two versions of the model, one 
elastic and the other inelastic. To get isospin violating couplings of the DM to nucleons, 
we give the DM an additional U(l) coupling to a vector Z' that kinetically mixes with 
the photon. 

The dark sector field content consists of two Weyl DM components xi,2 with charges 
ziz{gB,gz') under U(1)b xU(l)z', a real singlet a complex Higgs cp' with charge (0, 5^/), 
another <j) with charge {gsfi), and in the inelastic version of the model, a third (j) with 
charge {—2gB, —2gz') There are also "dark neutrinos" vi, 1^2 that carry no gauge quan- 
tum numbers. All of the Higgs fields need to get VEVs in order to give masses to the 
dark sector particles. The interactions are 

V = y^^xiX2 + yu'^i^ii^2 + M^U2J^2 + 9%\4>\^B^ +9'z^\4>'\^Z'^ 
+ {gsB^ + gz'Z'^){xWxi - A^^X2) + gB^$N + eep/p 
+ { ^{gBBi,^ gz''Z'^f\4>\^ ^\y^{4>x\X\^'\>*X2X2), inelastic model } (1) 

where p is the proton, and N is either kind of nucleon. Note that for the elastic version of 
the model, the last line is omitted, and expressions like X1X2 are shorthand for Xi"'^2X2- 
The lighter dark neutrino which we will call v' gets a small mass from the see-saw 
mechanism once ^ gets a VEV, due to the large bare mass My. The structure of the 
couplings of can be justified by a discrete symmetrjH xi ~^ ~Xi) ^1 ~^ ^ ~^ 
The light dark neutrino mass eigenstate i^' allows the DM to annihilate invisibly via 
XX — followed by ^ — )• z^V. For simplicity we have imposed a second discrete 
symmetry xi ^ X2t <P ^ (1^* so that there is only a single Yukawa coupling y. When 
$ and 4> get their VEVs, the DM mass eigenstates are x± = "^(xi =t X2) with masses 
M± = ± IX where = y^i^) and /i = y^{4>). 

The gauge bosons B and Z' do not mix with each other in the elastic version of 
the model, but the VEV of (j) causes them to do so in the inelastic version. We want 
to suppress such mixing to avoid having any gauge boson with significant couplings to 
both baryon and lepton number, since the constraints on such particles are quite severe. 
(These constraints also place stringent limits on isospin violating dark matter models in 
which a single vector couples to baryon number and mixes kinetically with the photon.) 
Ref. [27] shows that a vector that couples to B and has kinetic mixing of e = 10~^ 

^ This symmetry must be softly broken to avoid the cosmological domain waU problem, but such 
breaking in the potential for $ does not concern us in the following. 
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is relatively safe, but one with e = 10~^ is ruled out if its mass is less than 10 GeV. 
Moreover our determination of the correct amount of isospin violation is simpler if the 
mixing between B and Z' is negligible so that we can consider these fields to coincide 
with the mass eigenstates. We will thus assume that (0) is significantly smaller than 
{(j)) and {(j)'). In the mass eigenstate basis, the interaction of x with B becomes purely 
off-diagonal: qbX+^X-^ li.c, written in terms of Majorana-Dirac spinors. 

Fitting to CoGeNT/DAMA and relic density. We now consider the constraints 
arising from direct and indirect detection. (See table [T] for a summary of the results 
for constraints on all the model parameters.) Comparing the diagrams for x-it- ^-nd x-P 
scattering, the ratio of DM couplings to neutrons and protons is 

k = i + i^lf (2) 

In 9 b ^Z' 

We fix fp/ fn = —1.53, as needed to evade bounds from Xenon and to reconcile the 
DAMA/CoGeNT observations [TT]. The overall rates for DAMA/CoGeNT are repro- 
duced by matching the cross section for scattering on neutrons, an = dBl^n/ {^''^b) t° 
value 3 x 10~^^ cm^ |1HI15|. where the reduced nucleon mass is ^„ = 0.84 GeV if the DM 
mass is = 8 GeV. This gives niB/gB = 232.3 GeV and ml,/{gz'e) = -(79.9 GeV)^. 

As described in [15j, we can use the above relations between {mB,gB} and {mz',gz'} 
to eliminate the couplings in terms of the masses when computing the annihilation cross 
section into standard model particles by the processes XX ~^ BB, Z'Z' , BZ' — t- ffff or 
XX B,Z' ^ ff. For a given choice of e, this results in contours of (fXann^^) in the tub- 
mz> plane. By comparing the resulting value of ((Xannf ) to that required for the standard 
relic density, we see which ranges of the masses are compatible with annihilation that 
is subdominantly into visible particles, so that the relic density can be determined by 
invisible annihilations. The result is given in fig. [H which shows that for a given value 
of e, there exists a region where mzi and tub are both bounded from above such that 
annihilation into SM particles is sufficiently suppressed. For a given value of e, we find 
that the maximum allowed value of m^/ is given by 

mz' < \/e- 16.6 GeV (3) 

which corresponds to the horizontal part of the contours showing where {(Jaim 

v) < 0.1 • 

(<7ann^^)o in fig- [H This is the minimum reduction of annihilation into visible particles 
needed to satisfy the astrophysical constraints. On the other hand, < 4 GeV is the 
bound on regardless of e. 

To get the correct relic density, it is convenient to assume that < so that xX ~^ 
is allowedlll In this case we can determine a function of and m$ by demanding 
that the cross section for xX ~^ be 3 x 10~^^ cm^/s. The theoretical cross section 
is given by {a,nnv) = y^f{m^/M^)/{647rM^), where f{x) = (1 - x^/yil - x'^jl.f. 
For ?Ti$ < 0.9?Tii^, is in the range 0.076 to 0.11. This implies that (<I>) ~ 73 — 105 
GeV, since = y^{^) ^ 8 GeV is determined by the fit to CoGeNT/DAMA [llj. We 

"^If m# > M^, one must compute the annihilation cross section for xX '^v' . 
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Figure 1: For kinetic mixing e = 10~^, 10~^, 10~^, we show contours in the m^-m^' 
plane where the DM annihilation cross section into visible particles is equal to 
1 (solid), 0.1 (dashed) and 0.01 (dot-dashed) times the standard relic density 
value 1 pb-c. The allowed regions are below and to the left of the corresponding 
contour. 



find that the s-channel contribution to XX ~^ which is controlled by the cubic term 
in potential for <I>, is smaller than the dominant t-channel contribution by a factor of 
m|,/(16y^M^(<^)^) <^ 1 in the squared amplitude. Furthermore, annihilation to v'u' is 
suppressed because, as we show below, the Yukawa coupling y^, is much smaller than y^. 

New Higgs boson parameters. Next we turn to properties of the other new Higgs 
bosons and Yukawa couplings. For the inelastic version of the model, we require that 
gBgz'{4')'^ "ifi; '^l' to suppress the gauge boson mixing. If the gauge couplings are 
of similar size, this implies {(f)) <^ (0), {(j)'). Since the mass splitting between the DM 
states is given by 2/i = 2y^{(l)) = 10 keV, having small {(j)) helps to explain the smallness 
of this splitting. For example if {(f)) ~ 10 GeV, then y^ ~ 10~^. 

The other VEVs {(f)), ((/)'), can be be relatively large, if the couplings qb, g'z' are small. 
Using the constraints on {mB,gB) and {mz',gz') found above, we have 



{4>) = ^- 232 GeV, (0') = . 80 GeV (4) 

gn gz' 

In particular, by choosing to be smaller than gB, 4' can have a larger VEV than the 
SM Higgs. This helps the model to avoid being marginalized by recent very stringent 
constraints on the fourth generation of quarks needed for cancellation of the anomaly of 
gauged baryon number, since {(f>) can contribute to the mass to vector-like 4th generation 
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quarks [26]. We discuss this issue further in the final section. 

In the potential ([1]) we have omitted the renormalizable couplings of (p and ^ to the 
standard model Higgs fieldjl 

x^h\(i>\^h^+x^hm^h^ (5) 

which would give rise to (p-h and ^-h mixing, with the mixing angles 

f = — 2 2 ' ^ ~ — 2 2"' v"-* 

where v = {h) = 246 GeV. Such mixing gives rise to additional interactions between 
the DM and nucleons by exchange of the new Higgs particles. Because y,^ is so small, 
there is no strong constraint on 6, but $ exchange is potentially dangerous since it would 
perturb the amount of isospin violation away from the optimal value if its effect on x- 
nucleon scattering was comparable to that of the gauge boson exchange. Moreover the 
inelastic feature would be spoiled by strong $ exchange since the couplings of $ are not 
off-diagonal in the x mass eigenbasis. To avoid this problem we require that 

^ « 4 (7) 

where y„ = 0.3m„/f is the coupling of h to the nucleon (see for example [28]). However 
using the above values for gs/nT-B and y^., this is not very stringent; even for m$ = 4 
GeV it only constrains O <C 2.6, which far weaker than direct experimental constraints 
on the mixing of such a light Higgs boson, Q < 10~^ [29]. 

Constraints on dark neutrinos. Finally we come to the properties of the light dark 
neutrino, u'. It was last produced when the temperature of the universe was 0(m$) 
from the decays $ — t- i/V. If m$ ~ 1 GeV, this occurred before the QCD phase 
transition, and the relic density of z/' is suppressed by a factor of ~ 60, the reduction in 
the number of degrees of freedom in the thermal plasma. This means that cosmological 
upper bounds on the sum of neutrino masses, which we take to be 0.3 eV (see [30] for 
a review), are relaxed by a factor of ~ 22 (using the fact that ordinary neutrinos have 
a number density that is 4/11 that of photons), giving the bound niu' < 6.5 eV. This 
puts an upper limit on y^, 

y. < 3.5 X 10-^ j (8) 

Because of the dilution of the u' number density, there is no constraint from its contri- 
bution to the Hubble rate during BBN. 

A potentially serious constraint arises from the fact that i^' can be produced in su- 
pernovae by the nucleon collisions NN — )■ NNu'u' (via virtual <I> emission), causing 



^(j)' of course also has such couphngs but is unimportant to our discussion because it does not couple 
to the DM. 
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them to cool more quickly than observed in SN 1987A. A similar process in which ster- 
ile neutrinos are produced by vector current couplings to nuclei was considered in ref. 
[31]. (The vector coupling is expected to give similar results to our scalar coupling [32].) 
Adapting these results to the present case, we find that 

m$ ~ 100 TeV ^ ' 

Using our fiducial values m$ ~ 1 GeV and ~ 3.5 x 10^^, this translates into the 
constraint B < 2 x 10~^, which puts the mixing of <I> with h out of experimental reach 
[29] unless yy is taken to be smaller. 

For DM to annihilate invisibly, the decay channel — t- u'v' must dominate over 
decays into standard model channels like due to the mixing B, implying that 

y1 » (By^)^ (where y^ = 4.3 x 10~^ is the SM Yukawa coupling of /x). If B saturates 
its direct experimental upper bound of 10~^ and y^ saturates ([8]) then this criterion is 
satisfied even for My at the TeV scale. Heavier decay products can be accommodated 
by a modest increase in My. 

Interestingly, positive evidence for an additional species of dark neutrinos has come 
from recent CMB data [33j. If m$ is below the QCD scale, the decays V can occur 
sufficiently late for v' to have as large an abundance as the SM species. Smaller values 
of would naturally be correlated with lower v' masses since both are related to {^). 
These neutrinos must be sufficiently light to constitute an extra component of radiation 
before recombination, rather than being an extra component of the dark matter. 

Discussion. We have presented what we believe to be the simplest model of symmet- 
ric light dark matter that has the potential for explaining the tentative evidence from 
CoGeNT and DAMaIiI while robustly evading stringent astrophysical constraints that 
may soon exclude all such models unless the dominant DM annihilations are into invisi- 
ble particles. The isospin violating interactions between DM and nucleons are mediated 
by light vector bosons B of gauged baryon number and a Z' that has kinetic mixing e 
to the photon. 

The model is strongly constrained, and predicts a maximum value of mz' for each 
value of e ([3]) to sufficiently suppress the visible annihilation channels, as well as a 
minimum value of to suppress the gauge boson mixing angle. There are shown in 
fig. [2j The region of the m^/-e plane to the left of the solid line is interesting from the 
point of view of the proposed Heavy Photon Search (HPS) [M] and DarkLight [35l [36] 
experiments at Jefferson Laboratory. Fig. [2] shows that these two experiments, as well 
as the already-running APEX experiment [37] . are capable of probing a large fraction of 
the parameter space predicted by the model. We make the interesting observation that 
both B and Z' could be discovered by these experiments in the inelastic version of the 
model, where B has a small mixing 6b with Z' and thus couples to the electromagnetic 
current with strength eOs- 

®in general, the heaviest SM fermions that can be produced by "I> decay, keeping in mind that < 
M^'^S GeV 

''from footnote 1, it appears that that CRESST can only be made compatible with these if the require- 
ment of inelastic scattering is dropped 
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value 


constraint 


— e 


< 0.03 


fig.E] 


niB 


< 4 GeV 


fig.m 


niz' 


< • 16.6 GeV 


eq. © 


9B 


mB/(232 GeV) 


eq. ([2]), an 


gz' 


-e-im|,/(80 GeV)2 


eq. ([2]), C7„ 


(0) 




gauge boson mixing 


{4>) 


{gBrOB) ■ 232 GeV 


eq. (BD 


{4>') 


(\/bz'e/<7z')-80 GeV 


eq. © 


($) 


73 - 105 GeV 




Vx 


0.076 -0.11 


M^, relic abundance 




~ 10-^ 


X mass splitting 




< 3.5 X 10-^ 






~ TeV 





Table 1: Summary of allowed values of parameters of the model ([T]), and the correspond- 
ing constraints which determine them. 

Another experimental test of the model will come from the requirement of a fourth 
generation of quarks carrying exotic baryon number ±1 in order to cancel gauge anoma- 
lies of B |25[ ES] • Because the field cj) that spontaneously breaks B can have a larger 
VEV than the SM Higgs, it can increase the mass of vector-like exotic quarks without 
requiring them to have very large Yukawa couplings. Otherwise, such exotic quarks 
would be difficult to hide from current collider searches. The strongest constraint comes 
from the enhancement by a factor of 9 of Higgs boson production by gluon fusion [38] 
due to the extra generation. Tevatron and ATLAS have recently produced similar pre- 
liminary upper limits of m/^ < 124 GeV [39j and mh < 120 GeV [40j at 95% c.l. through 
the decay channel h — )• W~^W^ . Our model avoids these constraints to the extent that 
((/)) = (93/93) • 232 GeV can exceed the SM Higgs VEV. Barring a large hierarchy 
between these gauge couplings, it is unnatural to make (0) very large, so one may an- 
ticipate a softening of this constraint rather than completely evading it. In addition, 
the new Higgs boson (p can have large mixing with the SM Higgs, which would lead to 
further softening of the constraint. 

The model we have presented ties current hints for light dark matter to light vector 
boson searches, Higgs physics, fourth generation quarks, and the cosmology and astro- 
physics of sterile (dark) neutrinos. Even if this particular model should be ruled out by 
future data, the idea of invisibly annihilating dark matter may prove useful if its mass 
is at the 10 GeV scale. 

Acknowledgments. We thank C. Burgess, R. Essig, J. Lavalle, J. Kamenik, G. 
Raffelt, D. Spolyar, M. Trott, M. Williams, M. Wise, J. Zupan, for valuable discussions or 
correspondence. JC thanks the Aspen Center for Physics for a stimulating environment 
while this work was being done. 
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Figure 2: Solid (blue) line shows the maximum allowed value of m^/ (eq. ([3])) in the 
dark matter model for a given value of e, in the mz'-e'^ plane; dark regions 
are already excluded by various searches for light vector bosons, and light 
circumscribed regions are targeted by upcoming experiments. Background 
figure courtesy of R. Essig. 
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